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Abstract. We provide a data structure for maintaining an embedding of a graph on a surface (represented combi- 
natorially by a permutation of edges around each vertex) and computing generators of the fundamental group of 
the surface, in amortized time 0{logn + log g (log log g)^) per update on a surface of genus g; we can also test 
orientability of the surface in the same time, and maintain the minimum and maximum spanning tree of the graph in 
time 0{logn + log^ g) per update. Our data structure allows edge insertion and deletion as well as the dual opera- 
tions; these operations may implicitly change the genus of the embedding surface. We apply similar ideas to improve 
the constant factor in a separator theorem for low-genus graphs, and to find in linear time a tree-decomposition of 
low-genus low-diameter graphs. 

1 Introduction 

In this paper we introduce a tool, the tree-cotree decomposition, for constructing generators of fundamental 
groups of surfaces, and apply it to several algorithmic problems for graphs embedded on surfaces. 

The fundamental group of a surface measures the ability to transform one curve continuously into an- 
other. Two closed curves cq and ci on a surface are homotopic if there exists a continuous map from the 
cylinder [0, 1] x to the surface such that each cylinder boundary / x is mapped homeomorphically 
to curve c,. The elements of the fundamental group are equivalence classes of curves under homotopy. The 
identity element is the class of curves homotopic to the boundary of a disk, the inverse of a curve is formed 
by tracing the same curve backwards, and the concatenation of two curves forms their product in the group. 
For further exposition see any topology text, e.g. JT?]]. 

Fundamental groups are often infinite, but can be specified finitely by a system of generators and rela- 
tions. The generators of a fundamental group are a system of closed curves with the property that any other 
curve can be generated (up to homotopy) by concatenations of closed curves in the system. For instance, on 
a torus, two natural closed curves to use as generators are an equatoral circle gQ and a longitudinal circle 
g\; any other closed curve c on the torus is homotopic to a curve g^^^ that winds / times around the equa- 
tor followed by j turns around a longitude. The relations of a fundamental group describe nuU-homotopic 
combinations of generators: in the torus, the relation go^l^o '^F^ ~ ^ implies that the group is commutative. 

In addition to their topological interest, fundamental groups and their generators have other applications. 
Cutting a surface along its generator curves produces a planar surface with holes at the cuts, and Erickson 
and Har-Peled ^ have recently investigated the problem of finding short cutsets in nonplanar surfaces. As 
we will see, this cutting technique allows one to apply planar graph algorithms such as separator construc- 
tion [jl]| and tree decomposition to nonplanar graphs. We will also show how to use the generators to test 
orientability of surfaces and maintain minimum spanning trees of graphs embedded in surfaces. 

Subgroups of the fundamental group correspond to covering spaces, the most important of which is the 
universal cover. In the example of the torus, the universal cover is a plane, and the endpoints of curves ^q^^ 
lift to a planar lattice. If one chooses any pair of generators of the lattice, curves from the origin to these 
points map to a pair of generators for the fundamental group of the torus. Thus construction of generators is 
also closely related to lattice basis reduction, Euclid's algorithm, and continued fractions. 

1.1 New Results 

We show the following results. 



- We provide a simple dynamic graph data structure for maintaining a dynamic graph, embedded in a 
2-manifold, subject to updates that insert or delete edges, changing the ordering of edges around each 
vertex and thus implicitly changing the underlying 2-manifold. Our structure can also handle the dual 
operations of edge contraction and expansion, and maintains a set of generators for the fundamental 
group of the surface, as well as additional information such as the minimum spanning tree of the graph 
or the orientability of the surface, in amortized time 0{glogn) per update where g is the genus of the 
surface at the time of each update. 

- By combining our simple structure with ideas of separator based sparsification and recent poly- 
logarithmic algorithms for dynamic connectivity in general graphs [|^, 18 1 we further improve our 
time bounds, to 0(logn + log g (log log g)^) per update for generators and orientability, and 0(logn + 
(logg)^) per update for minimum spanning trees. 

- We improve by a factor of \/2 the constant factor in the best previous separator theorem for bounded 
genus graphs [|]]. 

- We provide an efficient algorithm for our previous nonconstructive result, that genus g graphs with 
diameter D have a tree-decomposition with ti"eewidth 0{gD) [^]. 



1.2 Related Work 

The previous work most close to ours is on construction of canonical schemata and cutsets. 

A canonical schema for a 2-manifold is a set of generators for its fundamental group, having a single 
prespecified relation ([{ g2ig2i+ig2i ^2(+l oriented surfaces, Y{gf for unoriented surfaces). The genera- 
tors must have a common basepoint, but unlike the ones in our construction they are not required to follow 



the edges of a given graph on the surface: they may pass across the interior of cells. Vegter and Yap | |19| ] 
showed that canonical schemata exist with total complexity 0{ng) where g is the surface genus; this bound 
is tight in the worst case, and they also showed that a canonical schema can be found in time 0{ng). The 
generators we find can have the same 0{gn) total complexity, but we can find an implicit representation of 
the generators in linear time. However, the generators we find are not necessarily in canonical form. 

Recently, Erickson and Har-Peled studied the problem of finding a cutset, that is, a set of edges 
the complement of which forms a topological disk on the given surfaces. Erickson and Har-Peled found 
algorithms for computing minimum-length cutsets in polynomial time for bounded-genus embeddings, but 
showed the problem to be NP-complete in the case of unbounded genus. In the applications to tree decom- 
position and separator theorems, we use our techniques to find cutsets with a guaranteed bound on the total 
length, but our cutsets may be far from the minimum possible length. 

Planar separator theorems have been long studied Jl4| ] and there has been much work on improved 



constants for such theorems Q16|]. More generally, for graphs embedded on genus g surfaces, a separator 



with 0{yjgn) vertices can be found in 0{n) time [[T]]. Specifically, Aleksandrov and Djidjev [|l|] showed 
that, for any e > 0, one can find a set of \/{l6g + 0{\/€))n vertices the removal of which partitions the 
graph into components of size at most en. Our techniques improve the constant factors in this result. We 
also apply similar ideas to the computation of tree-decompositions for low-diameter low-genus graphs; our 
previous work [§] proved the existence of such decompositions but did not provide an efficient algorithm 
for finding them. 

The main inspiration for this paper was previous work on dynamic connectivity and minimum spanning 
trees in embedded planar graphs which worked by decomposing the graph into two subgraphs, a min- 
imum spanning tree and a maximum spanning cotree. We use a similar decomposition here, however for 
surfaces of nonzero genus we end up with leftover edges that belong neither to the tree nor to the cotree. 
These extra edges are what we use to form the generators of the surface. We also combine the dynamic graph 
algorithms from that previous work, with a technique for maintaining small representative subgraphs [0] and 
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Fig. 1. Left: an embedding of 3 in the projective plane. The edges drawn crossing the dashed circle 
connect to the opposite side of the circle; there is one hexagonal cell in the center surrounded by three 
quadrilaterals that form a Mobius strip. Right: a gem representation of the embedding. The flags of the 
embedding are shown as triangles; the dashed edges represent connections of type Sy, the dotted edges 
represent connections of type Sg, and the soUd edges represent connections of type Sc- 



more recent methods of Holm et al. [12] to achieve time bounds that match the best of these dynamic graph 
algorithms. 

2 Representing an Embedding 

We consider maps: graphs embedded on two-dimensional manifolds, in such a way that the cells of the 
embedding are disks. We do not require that our graphs be simple. Such an embedding on an orientable 



manifold can be represented by specifying a clockwise ordering of the edges around each vertex [ |10| , Sed- 
ction 3.2.4]. However for our purposes it is convenient to use a somewhat more general representation, 
the graph-encoded map 13|, which also allows unorientable embeddings. We now briefly describe this 



representation. 

Suppose we are given a map, such as the embedding of ^3 3 shown on the left of Figure [I[ A flag of the 
map is defined to be the triple formed by a mutually adjacent vertex, edge, and cell of the map. For a flag F, 
let V{F) denote the vertex, E{F) denote the edge, and C{F) denote the cell. We define three operations that 
change one flag F into another: 

- Sv{F) is a flag formed by the edge and cell of F but with a different vertex, found at the other endpoint 
of E{F). 

- Se{F) is a flag formed by the vertex and cell of F but with a different edge, found by tracing around the 
boundary of C{F) from E{F) in the direction of V{F). 

- Sc{F) is a flag formed by the vertex and edge of F but with a different cell, the cell that is the next one 
clockwise or counterclockwise around V{F) from C{F) in the direction of E{F). 

Each of these operations produces a unique well-defined flag, and the neighboring relation between 
flags is symmetric: Sv{Sv{F)) = Se{Se{F)) = Sc{Sc{F)) = F. In a computer, one can represent these 
flags as objects in an object-oriented language, and store with each flag a pointer or reference to its three 
types of neighbor. Graphically, we view the flags as vertices in a graph, and the operations as unoriented 
edges, connecting each flag to its neighbors SviF), Se{F), and Sc{F), as shown on the right of Figure |l|. 
Although we see this gem as drawn in a projective plane, the gem representation of an embedding is purely 
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graph-theoretic and independent of any drawing of the gem. The edges of this graph should be viewed as 
having three colors, representing the type of operation generating the edge; in the figure these are shown as 
dashed, dotted, and solid, respectively. 

For each vertex of the map, with d adjacent edges, the resulting gem has a Id-cycle, with edges alternating 
between types Se and Sp- For each edge of the map, the gem has a quadrilateral with edges alternating 
between types Sy and Sq- And for each ^-gon cell of the map, the gem has a 2^-cycle with edges alternating 
between types Sy and Se- In this way, every map generates a 3-regular 3-edge-colored graph in which the 
Sy-Sc cycles all have four edges. Conversely, given any such graph, we can recover the original map by 
forming a ^-gon cell for each 2/:-cycle of edges of types Sy and Se and by gluing these cells together into a 
manifold whenever the corresponding edges of types Sy are part of a four-cycle of types Sy and Sc- 

Any map M. has a dual map A^* on the same surface, formed by creating a dual vertex/* within each 
face / of the primal map, and creating a dual edge e* for every primal edge e, so that if e is adjacent to two 
faces /i and/2, then e* connects /j* and/2 a path that crosses e once and crosses no other primal or dual 
edge. In the gem representation, the dual can be formed very simply, by reversing the roles of the gem edges 
of types Sy and Sc- 

In our dynamic graph data structures, we will augment the gem representation by balanced binary trees 
for the cycles corresponding to each vertex and cell of the map; in this way we can quickly look up which 
vertex's or cell's cycle contains a given edge of the gem, and whether two edges on the same cycle are in the 
same or opposite orientations around the cycle. 

3 Trees and Generators 

A spanning tree of a graph or map is just a tree formed by some subset of the edges of the graph that 
incorporates all the vertices of the graph. If C* is a spanning tree of the dual of a map, we call C = {e \ e* ^ 
C*} a spanning cotree of the map. If the edges of a map are given weights, the weight of a tree or cotree is 
defined to be the sum of the weights of its edges. 

Lemma 1. Let a map M be given, with distinct weights on each of its edges- Then the minimum weight 
spanning tree of M and the maximum weight spanning cotree of M are disjoint- 

Proof. Let e be an edge in the minimum spanning tree of . Then removing e from the minimum spanning 
tree results in a forest F with two connected components; let Fi be one of those components. Then the 
sequence of faces of Ai surrounding Fi forms a (possibly non-simple) cycle in A4* consisting of the duals 
of all edges connecting Fi to F2- Since e belongs to the minimum spanning tree, e* must be the shortest 
edge in this dual cycle, and so cannot belong to the maximum spanning tree of Ai*- □ 

In the special case of planar graphs, the minimum spanning tree and maximum spanning cotree form 
a partition of all the edges of the graph [^]. We define a tree-cotree partition of 7W to be a triple {T, C,X) 
where T is a spanning tree of Ai, C is a spanning cotree of Ai, and the three sets T, C, and X are disjoint 
and together include all edges of . In particular, if T is the minimum spanning tree and C is the maxi- 
mum spanning cotree, then {T, C,E{M.) \ {T U C)) is a tree-cotree decomposition. By assigning weights 
appropriately, we can use Lemma [I] to find a tree-cotree decomposition involving any given tree T or cotree 
C. More generally, if T and C* are forests such that T and C are disjoint, we can assign weights to make T 
become part of the minimum spanning tree and C become part of the maximum spanning cotree, extending 
T and C to a tree-cotree decomposition. 

We now show how these decompositions are connected to fundamental groups of surfaces. The funda- 
mental group of a space S is most commonly defined relative to some base point xq G S, but for suitable 
spaces (including the 2-manifolds considered here) it is independent of the choice of this base point. Define 
a loop to be a continuous function/ : [0, 1] 1— > 5 satisfying /(O) =/(l) = xq, and define two loops /q and/i 
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to be homotopy equivalent if there exists a continuous function/ : [0, 1]^ i-^ 5 such that/(x, /) = fi{x) and 
/(O, i) =/(!, /) = xq. Define the product of two loops /q and/j to be the loop 



/o(2x) forO <x < 1/2 
/i(2x-l) forl/2<x<l 



Then it can be shown that homotopy equivalence is an equivalence relation, and that the homotopy equiv- 
alence classes of loops form a group, the. fundamental group, under homotopy equivalence. In this group, 
inverses can be found by reversing the direction of a loop: if / is a loop, g{x) =/(l — x) is a loop inverse to 
/. A set of loops is said to generate the fundamental group if every other equivalence class of loops can be 
reached by products of these generators and their inverses. 

If we are given a rooted tree T and oriented edge e ^ T, with the base point xq of the fundamental loop 
at the root of T, we can define a loop loop(r, e) by following a path in T from xq to the head of e, then 
traversing edge e to its tail, and finally following a path in T from that tail back to xq. When considering 
these loops as generators of the fundamental group, the orientation of e is unimportant, since reversing the 
edge merely produces the inverse group element. 

Lemma 2. Let {T,C,X) be a tree-cotree decomposition of map M. Then the loops {loop(r, e) | e G X} 
generate the fundamental group of the surface on which M is embedded. 

Proof. Contract the edges in T into the root xq of T, while leaving the surface unchanged outside of a small 
neighborhood of T; this contraction does not change the fundamental group of the surface. In the contracted 
image, the loops from the statement of the lemma are each contracted into a single edge connecting xq to 
itself. Decompose the surface into a set consisting of a small neighborhood of the set of these contracted 
loops, and a set ^2 consisting of the faces of the map and the cotree edges in C. Then, by the Seifert - Van 



Kampen theorem []17[ Section 3.4] the fundamental group of the overall surface can be formed by combining 
the generators and relations of the fundamental groups of these two pieces, with additional relations describ- 
ing the way these two pieces fit together. S\ can be contracted to a graph with one vertex and \X\ self-loops, 
which has a fundamental group with each of these loops as generators and no relations. ^2 is topologically 
a disk (since it consists of a collection of disk faces glued together in the pattern of a tree) so it has a trivial 
fundamental group. The fundamental group of the overall surface, then, is generated by the loops described 
in the statement of the lemma, with a single relation formed by the concatenation of these loops around the 
boundary of ^2. □ 

Of course, for any tree T, the loops {loop(r, \ e ^ T} also generate the fundamental group, but 
in general this will form a much larger system of generators. The advantage of the system of generators 
provided by Lemma ^ is its small cardinality, proportional to the genus of the surface. 

The same tree-contraction argument used in Lemma § also shows that the edges in the loops {loop(r, e) \ 
e G X} form a cutset. 



4 Dynamic Generators 
4.1 Update Operations 

Following [^, we allow four types of update operation: insertion of a new edge between two vertices, 
deletion of an edge, dual edge insertion (expansion), and dual edge deletion (contraction). The effect of 
these updates on the gem representation is shown in Figure |2[ 

The location of an edge insertion can be specified by a pair of gem vertices p and q, shown as shaded 
triangles in the figure. We then subdivide the two gem edges Ssip) and SE{q) into paths of Se and Sq edges, 
preserving the connectivity of the Se-Sq cycles representing vertices in the gem, and create a new graph 
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Fig. 2. Effects of edge insertion and deletion and their duals on the gem representation. 

edge (represented in the gem by a 4-cycle of Sy and Sc edges) using the Sc edges of the subdivided paths, 
in such a way that the new gem edges SeIp) and Sgi^j) are adjacent to a common edge of type Sy- 

In the most typical case of an edge insertion (Fig. ^, upper right), the two edges SeIp) and SE{g) are 
distinct; each of these edges is subdivided into a three-edge path, and the two new vertices on each path are 
connected in pairs by new Sy edges. We distinguish three subcases of this operation, depending on which 
cells of the map contain the endpoints of the inserted edge: 

- In a cell-merging insertion, the two endpoints of the inserted edge are attached into distinct cells of 
the map; that is, the two gem edges Sgip) and Ssiq) belong to distinct Ss-Sy cycles. The effect of the 
insertion is to split each of these two cycles into paths, and to splice these two paths together into a 
single cycle via the newly created Sy edges. We are increasing the number of edges and decreasing the 
number of cells, so the surface's genus increases. Topologically, this can be viewed as gluing a handle 
between the two cells, and routing the new edge along this handle. 

- In a cell-splitting insertion, the two endpoints of the inserted edge belong to a single cell of the map, and 
if the Sg-Sy cycle representing this cell is consistently oriented, then SeIp) and Sgiq) both point towards 
or both point away from their endpoints p and q respectively. The effect of the insertion is to split the 
Sg-Sy cycle into two paths, and then reconnect the paths into two cycles, forming two cells connected 
across the new edge. We are increasing both the number of edges and the number of cells, so the genus 
and the topology of the surface are preserved. This was the only type of edge insertion allowed in [^]. 

- In a cell-twisting insertion, the two endpoints of the inserted edge belong to a single cell of the map, but 
SeIp) and Ssiq) are not consistently oriented. The effect of the insertion is to split the SgSy cycle into 
two paths, and then reconnect the paths into a single cycle, with the orientation of the two paths reversed 
from their previous situation. We are increasing the number of edges while preserving the number of 
cells, so the genus increases. Topologically, this can be viewed as gluing a crosscap into the cell, and 
routing the new edge through it. 

Each of these types of insertion has a complementary deletion, specified by the edge to be deleted, which 
we call a cell-splitting deletion, a cell-merging deletion, or an untwisting deletion respectively. In addition. 
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we allow the dual operations to insertions and deletions, as shown on the upper right of the figure; dual 
deletion is more familiar as the operation of edge contraction. The dual of a cell-merging insertion merges 
two vertices of the graph while increasing the genus, the dual of a cell-splitting insertion splits a vertex 
into two adjacent vertices, and the dual of a cell-twisting insertion changes the order of the edges and cells 
around a vertex. We do not describe these in detail because our data structure will be self-dual, so that each 
dual operation is performed exactly as the primal operation substituting Sy for Sc and vice versa. 

Finally, it is possible to specify an insertion in which Ssip) and Ssiq) are not distinct. If p q while 
^sip) = Ssiq) (Fig. |2[ lower left), we subdivide Se into a path of five edges, and connect the four new gem 
vertices into an SySc cycle in such a way that adjacent vertices of the cycle are connected to p and q. This 
can be viewed as a special case of a cell-splitting insertion in which the two endpoints of the inserted edge 
are placed next to each other at the same vertex, creating a one-sided cell contained within the previously 
existing cell. The dual of this operation adds a new degree-one vertex to the graph. 

lfp = q, we perform an edge insertion by again subdividing Se into a path of five edges, and connecting 
the four new gem vertices into an SySc cycle in such a way that opposite vertices of the cycle are connected 
to p and the other endpoint of SeIp)- Edge insertion with p = q and dual edge insertion with p = q produce 
the same result, shown in the lower right of Fig. This can be viewed as a special case of a cell-twisting 
insertion in which the two endpoints of the inserted edge are placed next to each other at the same vertex. 



4.2 Simple Dynamic Generators 

We now describe a simple scheme for maintaining generators and minimum spanning trees in dynamic maps, 
based on the dynamic planar graph algorithms of Eppstein et al. [0]. We will later show how to improve the 
time bounds for this scheme at the expense of some additional complexity. 

Theorem 1. Let M be a map, with distinct edge weights, changing according to the update operations 
described in the previous section as well as by changes to the weights of the edges. Then in time 0{g\ogn) 
per update we can maintain a tree-cotree decomposition (T, C,X), where T is the minimum spanning tree, 
C is the maximum spanning cotree, and g is the genus of the map at the time of the update operation. 



Proof. We represent T by Sleator and Tarjan's dynamic tree data structure []15[], modified as in [J7|] to support 
edge contraction and expansion operations as well as the more usual edge insertions and deletions. We use a 
similar data structure to represent the dual tree C* . With this data structure, we can find in logarithmic time 
the longest edge on any path of T, or the shortest edge on any path of C* . We can also test in logarithmic time 
whether a pair of edges {e, e'} forms a swap; that is, whether T/\{e, e'} (where A denotes the set-theoretic 
symmetric difference) is another spanning tree of T. We also maintain the tree data structures discussed in 
Section ^ for the adjacency lists of each vertex and face of ^A, which allow us to classify the type of an 
insertion or deletion in logarithmic time. 

To increase the weight of an edge e in T U X, we first test whether the increased weight causes the edge 
to move to the maximum cotree C, by finding the shortest edge e' on the dual path connecting the endpoints 
of e*, and testing whether the dual swap {e, e'} increases the weight of C. If so, we add e to C and move e' 
to X. We then search X for the minimum weight edge e" forming a swap {e, e"} in T, and testing whether 
this swap causes an improvement in the weight of T. If it does, we add e" to T and move e to X. Lemma |I] 
shows that, after these steps, we must again have a tree-cotree partition. To decrease the weight of an edge e 
in C U X, we perform a similar sequence of steps on the dual tree-cotree-partition (C*, T*,X*). 

We now describe the effect on our data structure of the various insertion and deletion operations. The 
edge expansions and contractions are performed by very similar sequences of steps in the dual map. 

- To perform a cell-merging insertion of an edge e, we find the lightest edge e' on the dual path in C* 
connecting the merged cells, and place both e and e' in X. We then perform swaps similar to those for a 
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weight change to determine whether e should move into T or C. To perform a cell-splitting deletion of 
e, we search X for the heaviest edge e' connecting the two components of C* formed by the spht, and 
move e' into C. 

- To perform a cell-splitting insertion of an edge e, we add e to X, search X for the heaviest edge e' 
connecting the two components of C* formed by the split, and perform swaps similar- to those for a 
weight change to determine whether to move e into T or C. To perform a cell-merging deletion of e, 
we increase the weight of e until it belongs to C, then remove it from the graph and merge its two dual 
endpoints. 

- To perform a cell-twisting insertion of an edge e, we simply add e to X, search for the best swap {e, e'} 
with e' G r, and perform this swap if it improves the weight of T . The inserted edge e can not belong to 
C since e* connects a dual vertex to itself. To perform an untwisting deletion, we test whether e belongs 
to r, and if so swap it with the best replacement e' ^X before removing it. 

Each of these operations maintains the desired partition via a sequence of 0{g) dynamic tree queries 
and updates, so the total time per change to is 0(g log n). □ 

If we wish to maintain a tree-cotree decomposition of a graph without weights, we can choose weights 
for the edges arbitrarily. 



4.3 Improved Dynamic Generators 

The bottleneck of Theorem |l] is the search for swaps in graphs T U X or C* U X* ; performing this search by 
testing each edge in X* takes time 0{g log n). We can improve this by using the method of separator based 
sparsification to replace T and C* by contracted trees T' and C' that match the 0{g) size of X, and then 
applying a general dynamic graph connectivity algorithm [|2, 18] to find the swaps in the smaller graphs 
r'UXor C'UX*. 

Specifically, we form T from T by repeated edge contractions, forming a sequence of trees Tq = 
T ,Ti, . . .Tk = T'. If Ti contains a vertex with degree one that is not an endpoint of an edge in X, we 
delete that vertex and its adjacent edge to form r,+i. Otherwise, if Tj contains a vertex with degree two that 
is not an endpoint of an edge in X, let e and e' be the two edges incident to that vertex; we form r,+i by 
contracting whichever of e and e' has the smaller weight. The construction of C' from C* is essentially the 
same except that we contract the larger weight of two edges incident to any degree two vertex. 

Lemma 3. The trees T' and C described above have 0(1X1) edges, and do not depend on the order in 
which the contractions are performed. 

Proof. Consider the set of paths in T connecting endpoints of X. The union of these paths forms a subtree T" 
of T. Then an alternative description of T' is that it consists of one edge for each path of degree-two vertices 
in T", with the weight of this edge equal to the maximum weight of an edge in the path. Specifically, each 
deletion of a degree one vertex removes an edge that is not in T", and each contraction of a degree two 
vertex reduces the length of one of the paths in T" without changing the maximum weight of an edge on the 
path. This proves the assertion that T' is independent of the contraction order; it has 0(|X|) edges because 
it has at most |X| leaves and no internal degree-two vertices. □ 



Lemma 4. Trees T' and C' can be updated after any change to T or C* in time 0{logn) per change. Each 
change to T, C, or X causes 0{\) edges to change in T' and C'. 



Proof. We use the following data structures for T: First, we use a dynamic tree data structure that can 



find the maximum weight edge on any path of T, and that can determine for any three query vertices of T 
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which vertex forms the branching point between the paths connecting the three vertex. Second, we maintain 
an Euler tour data structure [ [III ] that can find, for a query vertex in T, the one or two nearest vertices that 
belong to T' . Whenever a change to X adds a vertex v to the set of endpoints in X, we must also add v 
to T'; we do this by querying the Euler tour data structure to find the edge or vertex of T' at which the 
branch to v should be added, querying the dynamic tree data stiucture to find the coiTcct branching point on 
that edge, giving the overall graph structure of the new tree T' , and finally querying the dynamic tree data 
structure again to find the weights of the changed edges in T' . To remove a vertex from the endpoints in X, 
we simply contract one or two edges as necessary in T' . To perform a swap in T, we use the data structures 
as described above to make all swap endpoints part of T' , perform the swap, and then contract edges to 
remove unnecessary vertices from the changed T'. The updates to C' are similar. □ 

Lemma 5. Let e be an edge of T that was not contracted informing T', and let e' be an edge ofX. Then 
{e, e'} is a swap in T if and only if it is a swap in T'. 

Lemma 6. Let e' be an edge in X, and let e be the heaviest edge forming a swap {e,e'} in T. Then e is not 
contracted informing T'. 



Theorem 2. Let Ai be a map, with distinct edge weights, changing according to the update operations 
described in the previous section as well as by changes to the weights of the edges. Then in amortized 
time 0{logn + (logg)^) per update we can maintain a tree-cotree decomposition {T, C,X), where T is the 
minimum spanning tree, C is the maximum spanning cotree, and g is the genus of the map at the time of the 
update operation. 

Proof. The overall algorithm is the same as the one in Theorem |l], but we use the data structures described 
above to maintain trees T' and C' , and use the dynamic minimum spanning tree data structure of Holm et 



al [12| to find the edge e' G X forming the best swap for any edge e G T or e* G C*. The time bound follows 
from plugging the 0(1X1) bound on the number of edges in T' U X and C' U X* into the time bound for 
the minimum spanning tree data structure, and using this bound to replace the 0{g\ogn) time for finding 
replacement edges in Theorem |l|. □ 

The final improvement to our bounds comes from the observation that, if we only wish to maintain a 
tree-cotree decomposition for a graph without weights, then we can replace the dynamic minimum spanning 
tree data structures for T' U X and C' U X with any dynamic connectivity data structure that maintains a 
spanning tree of these two dynamic graphs and changes the tree by a single swap per update to the graphs. 



The connectivity structure of Holm et al [12], as improved by Thorup [18], is a suitable replacement 



Theorem 3. Let Ai be a map, with distinct edge weights, changing according to the update operations 
described in the previous section as well as by changes to the weights of the edges. Then in amortized time 
0(log« + logg(loglogg)-^) per update we can maintain a tree-cotree decomposition {T, C,X), where g is 
the genus of the map at the time of the update operation. 



4.4 Orientability 

The results above are nearly sufficient to determine the surface underlying the dynamic map Ai, since it 
is known that 2-dimensional closed surfaces are classified by two quantities: their genus (which we can 
calculate from |X|) and the orientability of the surface. We can compute the orientability from the following 
simple result, which follows from the fact that the loops generated by edges in X form a cutset: 
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Lemma 7. Let M have a tree-cotree decomposition (T, C,X). Then The surface on which Ai is embedded 
is unorientable if and only if every sufficiently small neighborhood of at least one of the loops {loop(r, e) \ 
e G X} is homeomorphic to a Mobius strip. 

Theorem 4. Let M be a map, with distinct edge weights, changing according to the update operations 
described in the previous section as well as by changes to the weights of the edges. Then in amortized time 
©(log n + log g(log log g)^) per update we can determine after each update whether the surface underlying 
the map is orientable. 

Proof. For any path in A^, we can form a sufficiently small neighborhood of the path by gluing together 
a collection of disks, where each disk has as its boundary the gem cycle representing a vertex or edge in 
the path. We choose an arbitrary clockwise orientation for the gem cycle representing each vertex in Ai, 
and augment the Sleator-Tarjan dynamic tree data structure used to maintain the tree T in the tree-cotree 
decomposition (r, C,X), by storing an extra bit per tree edge that is zero when the edge's preserves the 
chosen orientations of its two endpoints and one when it reverses these orientations; we can use the dynamic 
tree data structure to determine the parity of the sum of these bits along any path in T. After each change 
to the tree-cotree decomposition, we use this data structure to determine whether any e £ X forms a loop 
with odd total parity (including the parity of e itself). If any odd loop is found, the map is unorientable; if 
all loops are even, the map is orientable. □ 

5 Separator Theorems 

Since Lipton and Tarjan's discovery of the planar separator theorem [[l4|], separator theorems have become 
widely used in algorithms, and there has been much work on improved constants for these theorems. For- 
mally, we define an e-separator of an «-vertex graph G to be a subset S of the vertices of G such that each 
connected component of G\S has at most en vertices. For graphs embedded on a genus-g orientable surface, 
Aleksandrov and Djidjev showed that there always exists an e-separator with l^l < \/{l6g + 0{l/e))n. 
We combine our tree-cotree decomposition idea with the Lipton-Tarjan method of partitioning the breadth 
first search tree into levels (also used by Aleksandrov and Djidjev) to improve this result and extend it to 
unorientable surfaces. 

By Euler's formula, if {T, C,X) is a tree-cotree decomposition of a map on an orientable surface, then 
\X\ = 2g where g is the genus of the surface on which Ai is embedded. If the surface is unorientable, 
then \X\ = g. Because of this confusion between orientable and unorientable surfaces, we prefer to use 
the Euler characteristic = 1^1 ~ 2 in stating our results. In terms of x, Aleksandrov and Djidjev's 

result is that there exists a separator with l^l < ^/ (8x + 0{l/e))n. We improve this by a factor of \/2, to 
\S\ < V{4x + 0{\/e))n. 

Let r be a breadth first search tree of the given graph, and form a tree-cotree decomposition {T,C,X). 
For any pair of integers < / < ^, define vertex sets Lf j^ and 5, ,t as follows: L,- consists of all vertices v 
whose depth d{v) in T satisfies d{v) = i (mod k). 5,^ is formed by adding to L,- ^, a path Pj^ki^) for each 
endpoint v of an edge in X, where Pi^k{v) C T connects v to the first member of L, ,t on the path from v to 
the root of T. 

Lemma 8. For any i and k, G\ Sj^^ is planar 

Proof. Extend the depth function on vertices to a continuous function A : 7W M that maps each vertex 
to its depth, maps each edge of G one-to-one onto the interval between the depths of its endpoints, and 
maps each cell of M. in such a way that the map has no local minima or maxima in the interior of the 
cell. Partition the surface underlying the map Ai into submanifolds with boundary, by removing all points 
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X for which A(x) = / mod k. Then each submanifold is formed from M. by cutting certain curves at depth 
X{x) = i + ak (for some integer a) and certain other curves at depth A(x) = / + (a + l)k. If we cap off these 
cut curves with disks, we form a manifold Ma without boundary. Within each of the added disks at depth 
\{x) = i + ak, add a vertex, connected by edges to all the vertices on the boundary of the disk; the result is 
a graph Ga embedded mM.a, such that each edge and vertex in G \ 5/ yt belongs to exactly one of the graphs 
Ga- If we let Ta be formed from the union of T n Ga with the newly added edges in Ga, then Ta is a forest 
disjoint from the coforest C H Ga, so we can extend these two subgraphs to a tree-cotree decomposition 
(r^, C'a,X'^) such that X'^ C XnGa- Therefore, the paths Pi^k{v) for v an endpoint of an edge inXn G^ form 
a cutset for 7Wa> and 5, is a cutset for all of A4. □ 

Lemma 9. For any k, and any map embedded on a surface with Euler characteristic x> 

E \Si,k\<n + {x + 2)k{k-\). 
0<i<k 

Proof. Each vertex of the map belongs to exactly one graph Ljj^. In addition, each of the 2(x + 2) endpoints 
V of edges in X contributes J2 \Pi,ki'^)\ = J2Q<i<k ' vertices to the total. Some vertices may be in more than 
one path Pi^ki^) but this only reduces the total. □ 

Theorem 5. Any n-vertex graph G embedded on a surface with Euler characteristic x has an e-separator 
with \S\< ^(4x + 0{\/e))n vertices. 

Proof. By choosing / minimizing \Si^k\, we achieve \Si^k\ < n/k + x^ + 0{^)- The result follows by choosing 
^ = vWX' applying a planar separator theorem to the planar graph G \ Si^k, and forming the union of the 
planar separator with 5,- yt- □ 

If a map representing a characteristic-x embedding of G is given as input, the method described in the 
proof of the theorem can be implemented to run in time linear in the size of the map, and produces a cutset 
of size at most \/(4x + 0(1))?!. The time bound for the overall algorithm, and the constant factor in the 
0{n/e) term in the separator size, depend on the choice of planar separator algorithm used to complete 
the separator construction. 

We note that not all loops need always be incorporated in Si^k in order to make G \ 5, ,t planar. For 
instance, if G is embedded on a torus, it can be cut into a planar graph in the form of a cylinder by removing 
only one of the two loops. In addition, if we choose Si^k in such a way that some of the subgraphs Ta and 
C n Gq described in the proof of Lemma ^ are disconnected, then additional edges are removed from X in 
forming the tree-cotree decompositions of J\Aa and again we can reduce the number of loops. Can we use 
these observations to further reduce the constant factor in our separator theorem? 

6 Low-Diameter Tree Decomposition 

In previous worklf] we showed that in a minor-closed graph family, the treewidth of a graph with diameter 
D can be bounded by a function of D, if and only if the family excludes some apex graph. For instance, for 
planar graphs (which exclude the apex graph K^) the treewidth is at most 3D — 0(1) [^. Such bounded 
local treewidth results can be used to find efficient solutions to many important algorithmic problems on 
these graphs 

In particular, any graph embedded on a genus g = 0(1) surface has bounded local treewidth, since the 
family of such graphs is minor-closed and excludes an apex graph. However, the treewidth bound from this 
general result grows rapidly with g, and in this special case we were able to prove a tighter bound: 



II 



Theorem 6 (Eppstein [§]). Let G have genus g and diameter D. Then G has treewidth 0{gD). 

The proof relies on showing the existence of a cutset of size 0{gD), and forming a tree-decomposition 
by adding this cutset to each node in a tree-decomposition of the remaining planar graph. However, one step 
in the proof is to find the minimum cardinality cutset, so the proof is nonconstructive, and Erickson and 
Sar-Peled [j8|] later showed that in fact finding minimum cutsets is NP-complete. By replacing this step with 
methods based on our tree-cotree decomposition, we can turn this nonconstructive proof into an efficient 
algorithm: 

Theorem 7. Let G have genus g and diameter D, and suppose that we are given as input an embedding M 
of G on a surface of genus g. Then in time 0{gDn) we can find a tree-decomposition of G with treewidth 
0{gD). 

Proof. We form a tree-decomposition (T, C, X) by letting 7 be a breadth-first spanning tree of G, and letting 
C be any spanning cotree disjoint from T. We then let S be the set of vertices incident to the edges in loops 
{loop(r, e) I e G X}. Each loop has at most 2D edges, so \S\ = 0{gD), and G \ 5 is planar, so it has a tree- 
decomposition of treewidth 3D — 0(1) that can be found in time 0{Dn) We form a tree-decomposition 
of G by adding S to each node in the tree-decomposition of G\S. □ 

The 0{gDn) time bound is optimal in the sense that it is the size of an explicit representation of the 
resulting tree-decomposition, but if we represent the sets of vertices at each tree-decomposition node im- 
plicitly, as a union of a constant number of paths in the two trees T and B, the time can be reduced to linear 
in the size of the input map M.. 



7 Conclusions 



We have shown how the tree-cotree decomposition can be used in several applications for graphs embedded 
on surfaces: maintaining dynamic properties of the graph, computing separators, and constructing tree- 
decompositions of small treewidth. There are many other important planar graph algorithms that can be 
generalized to graphs embedded on surfaces, and it would be interesting to see which others of these can be 
improved by the use of our tree-cotree decomposition methods. 
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